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ABSTRACT 

We present Hubble Space Telescope UV spectra of the 4.6 h period double white dwarf 
SDSS J125733.63+542850.5. Combined with Sloan Digital Sky Survey optical data, 
these reveal that the massive white dwarf (secondary) has an effective temperature 
T 2 = 13030 ± 70 ± 150 K and a surface gravity log <72 = 8.73 ± 0.05 ± 0.05 (statistical 
and systematic uncertainties respectively), leading to a mass of M 2 = 1.06 M 0 . The 
temperature of the extremely low-mass white dwarf (primary) is substantially lower 
at Ti = 6400 ± 37 ± 50 K, while its surface gravity is poorly constrained by the data. 
The relative flux contribution of the two white dwarfs across the spectrum provides 
a radius ratio of R 1 /R 2 — 4.2, which, together with evolutionary models, allows us 
to calculate the cooling ages. The secondary massive white dwarf has a cooling age of 
~ 1 Gyr, while that of the primary low-mass white dwarf is likely to be much longer, 
possibly > 5 Gyrs, depending on its mass and the strength of chemical diffusion. These 
results unexpectedly suggest that the low-mass white dwarf formed long before the 
massive white dwarf, a puzzling discovery which poses a paradox for binary evolution. 

Key words: stars: individual: SDSS J125733.63+542850.5 - white dwarfs - binaries: 
close 


1 INTRODUCTION 


Double white dwa rf binaries are common end products o f 
binary evolution llMarsh et al.l 1 19951 : iToonen et al.l l2014h . 
Those with separations small enough to have experienced 
one or two common envelope phases are particularly in¬ 
teresting, as t hey are thought to be progenitors of super- 
nov ae Type la (llben fe Tutukovll 19841 : 1 Webbinkl 19841) , Type 


.la (iBildsten et al. 20071). R CrB stars (Webbink 1984) and 
AM CVn systems (iBreedt et al.ll2012l : iKilic et~ l2014bl ) . In 
addition, mergers o f Galactic double white dwarfs occur rel¬ 
atively frequently (iBadenes fe Mao3l2012l ). and constitute 
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the main source of the background gr avitational wave signa l 
at frequencies dete ctable from space (iNelemans et al.ll200ll : 
iHermes et al.ll2012l) . 

There is an important relation between the initial 
mass of a main-sequence star and the final mass of the 
white dwarf that will form the remnant of that star 
(jWeidemannll2Qo'ol ). This initial-final mass relation predicts 
that extremely low-mass (ELM) white dwarfs, typically with 
masses M w< j < 0.3 Mq, cannot yet form as a natural prod¬ 
uct of stellar evolution because the main-sequence lifetime 
of their low-mass progenitors is longer than the present age 
of our Galaxy. However, ELM white dwarfs can be formed 
in binary systems in which the separation is close enough 
for the two stars to interact significantly before the ELM 
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progenitor has evolved off the main-sequence (mass trans¬ 
fer via Case A or early Case B Roche-lobe overflow). The 
binary companion causes the evolution of the ELM pro¬ 
genitor to be truncated before ignition of helium, and af¬ 
ter ejection of the envelope the helium core is exposed 
as the ELM white dwarf. Typically, ELM white dwarfs 
have surface gravities log g < 7, as well as re lativel y mas¬ 
sive hy drogen envelopes (~ ICC 3 ICE 2 Mm: llstrate et aD 
l2014bl b New dedicated searches such as the ELM Survey 
have signif icantly increased the known population in re- 
cent years (I Brown et alJl201tl 12011 l20ll iKilic et al.lboi ll . 
120121 ). The majority of ELM white dwar fs are companions to 
other white dwarfs llKaplan et alJ 2014 , this paper) or mil¬ 
lisecond pulsars ( s ee for example vanKerkwiik et al. l l 19961: 
iBassa et al.l 120061 : lAntoniadis et al.1120131') a nd a few h ave 
been found in hierarchical triple systems dRansom et al.l 
120141 : IKilic et a l.||2014a| , 20151) or orbiting A- or F-type main- 


sequence stars ( Maxted et alj|20l3 : iBreton et al .iLOil i . The 
subject of this paper, SDSS J1257+5428, is a binary that 
likely belongs to the first of these classes, but, as we shall 
see, how it evolved into the system we see today is a mystery. 


1.1 Introduction to SDSS J1257+5428 

The double white dwarf binary SDSS J1257+5428 (full 
name: SDSS J125733.63+542850.5) was first discovered 
when the available Sloan Digital S ky Survey (SDSS; 
lEisenstein et alj|2006l : lYork et al Jl2()()0l ) subspectra were ex¬ 
amined for radial velocity variations as part of the Sloan 
White dwArf Radi al velocity Mining Survey (SWARMS; 
iBadenes et al.ll2009l l . Follow-up spectroscopy revealed radial 
velocity variations with a semi-amplitude of 323 km s~\ 
which were interpreted to come from a 0.9 Mg white dwarf. 
Combined with the orbital period of 4.6 h and the absence 
of additional spectral features, this suggested that the most 
likely companion wo uld be a neutron star or a black hole 
(iBadenes et al.ll2009l j. 

Follow-up B and R band spectroscopy revealed two dis¬ 
tinct components in the spectra, although th e Balmer ab¬ 
sorpt i on lines only showed a single sh arp core dMarsh et al.l 
l201ll ; iKulkarni fc van Kerkwiikll2010l j . These deep, radial- 
velocity variable Balmer lines in fact originate in a cool, 
ELM white dwarf, which we hereafter refer to as the primary 
(because it dominates the flux at visual wave lengths, and 
follow ing IKulkarni fc van Kerkwiikl l20ld and iMarsh et al.l 
1201 ll j. The secondary is another white dwarf, which is hot¬ 
ter and significantly more massive, causing it to have very 
broad absorption lines. In addition, it is likely rotating fast, 
causing its line cores to be smeared out. Due to the shal¬ 
low nature of these lines, and the absence of sharp cores, it 
was not possible to detect a radial velocity variation of the 
massive white dwarf. 

The combination of these two white dwarfs in the same 
binary system is very interesting. The primary component 
is of much lower mass, and therefore has a much larger sur¬ 
face area than the secondary component. This causes the 
cooler primary to dominate the flux at wavelengths A > 
4000 A. At shorter wavelengths the secondary white dwarf 
starts dominating due to its higher temperature. Note that 
the fact that the higher mass white dwarf is hotter is con¬ 
trary to expectation since it presumably formed much earlier 
than the low-mass white dwarf. At the time of the studies 


bv IKulkarni fc van KerkwifO (l20ld l and lMarsh et al.l (1201 llj 
there were only a limited number of low-mass white dwarf 
models available, leaving it unclear whether or not the cool, 
low-mass white dwarf could have overtaken the secondary 
white dwarf on the cooling track. To securely measure the 
secondary’s temperature, we have obtained Hubble Space 
Telescope far-ultraviolet spectra. These new measurements 
of the hot white dwarf are presented in this paper and com¬ 
bined with recent b i nary models for ELM hel ium white 
dwarfs (llstrate et alJ l2014bl : lAlthaus et al.l l2013l l to study 
this binary’s evolutionary history further. 


2 OBSERVATIONAL DATA 

2.1 The Hubble Space Telescope data 

SDSS J1257+5428 was observed with the Hubble Space Tele¬ 
scope ( HST ) in Cycle 18, with program ID 12207. Part of 
the observations were done with the Cosmic Origins Spec¬ 
trograph (COS) on 2011 May 9, with the G140L grating 
and a central wavelength of A ce n = 1280 A. The total ex¬ 
posure time of these data is 146 minutes. The double white 
dwarf was also observed with the Space Telescope Imaging 
Spectrograph (STIS), on 2011 Oct 22. For these observa¬ 
tions, totalling 95 minutes, the G230L grating was used at 
Acen = 2376 A. The raw data were processed by the standard 
pipeline at the Space Telescope Science Institute. 

In the following analysis we exclude parts of the HST 
spectra that are contaminated by geocoronal Ol (1304 A) 
emission. In addition, for the COS and STIS data, we 
have excluded data at wavelengths A > 1700 A, and A < 
1650 A respectively, where the signal-to-noise ratio is very 
low. The measured flux is consistent with the Swift Ultra- 
Violet / Optical Telescope data presented in IMarsh et al.l 
(| 2011 |) . 

2.2 Parallax observations 

We used the 2.4m Hiltner telescope at the MDM Observa¬ 
tory on Kitt Peak on 19 observing runs between 2010 Jan¬ 
uary and 2014 June. The astrometric solution includes 128 
exposures, all taken in the /-band. Observations, reductions, 
an d analysis fol l owed proce dures similar to those des cribed 
in IThorstensenl (l2003l j and iThorstensen et all (l2008lf . The 
parallax of SDSS 1257+5428 relative to the reference stars 
was 8.3 mas, with a formal fitting error of only 0.8 mas, al¬ 
though we judged the external error to be 1.3 mas from the 
scatter of the reference stars. The colours and magnitudes 
of the reference stars yield a 1.6 mas correction due to the 
finite distance of the stars forming the reference frame, so 
our absolute parallax estimate is 9.9+1.3 mas, which on face 
value gives a distance to SDSS J1257+5428 of ~ 101 + 15 pc. 
The proper motion relative to the reference frame is mod- 
est, [/xx,/xvl = [—4 5,+91 mas yr _1 ; the PPMXL catalogue 
dRoeser et al . I Soil gives 1—4 1 .0,+1 1.81 mas yr \ in very 
good agreement. IThorstensenl ll2003l l describes a Bayesian 
procedure for estimating a distance by combining paral¬ 
lax information with proper motion (interpreted using an 
assumed space-velocity distribution) and with photometric 
distances. For the present case, we used only the proper- 
motion constraint to avoid tautology. The small proper mo- 
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tion combines with the Lutz-Kelker correction to give an 
estimated distance of II 2 IJ 5 pc, slightly larger than the in¬ 
verse of the parallax. Assuming a thick-disk velocity distri¬ 
bution increases this by another ~ 5 pc. 


2.3 ULTRASPEC photometry 

On the nights of March 2 and March 3, 2015, we observed 
SDSS J1257+5428 wit h the high -speed photometric cam¬ 
era ULTRASPEC dDhillon et al.l[20141 '). which is mounted 
on the 2.4 metre Thai National Telescope located on Doi 
Inthanon, Thailand. In total, we obtained 240 minutes of g' 
band dat a. The data were re duced using the ULTRACAM 
pipeline dDhillon et alj 12007( 1 . with which we debiased and 
flatfielded the data and performed relative aperture pho¬ 
tometry using a nearby bright star to minimise the effects 
of atmospheric variations in the light curves. 


4 RESULTS 

For each of the free parameters, the mean value and la un¬ 
certainty of the converged MCMC chain are listed in TableQ] 
column 2. Note that the quoted uncertainties are purely 
statistical. They do not include any systematic uncertain¬ 
ties that may be present, and are therefore underestimates 
of the true uncertainties. For our best model, the reduced 
xf = xf — 0.5. However, scaling the errorbars on the data 
such that Xv — 1 would only decrease the statistical un¬ 
certainties further, and we refrain from doing so. The re¬ 
sults from our MCMC are shown in Fig. [l] projected on 
the various 2 -dimensional parameter planes, as well as in 1 - 
dimensional histograms. The best model, together with the 
HST and SDSS data, is shown in Fig. [2] and fits the data 
well at all wavelengths. The underpredictions of the model 
with respect to the u and g SDSS fluxes (shown in the bot¬ 
tom panel) are less than 3cr of the SDSS flux. Given that 
SDSS uncertainties do not include systematic uncertainties, 
we do not think this difference is significant. 


3 FITTING SPECTRA USING A MARKOV - 

CHAIN MONTE CARLO APPROACH 

We fit both the HST COS and STIS spectra as well as 
the SDSS ugriz fluxes with a Markov Chain Monte Carlo 
(MCMC) analysis, using the affine-i nvariant ensemble sam¬ 
pler i n the python package EMCEE dForeman-Mackev et al.l 
120131 1. To obtain the SDSS fluxes, we use the PSF magni¬ 
tudes, which we correct for the offset between the SDSS and 
AB magnitude systems using (-0.04, 0, 0, 0, 0.02) for the 
ugriz measurements respectiveljQ. 

We fit t he data w ith a sum of two white dwarf model 
spectra from lKoesterl (|20ld l. which employ a mixing length 
ML2/a = 0.8, and list the Eddington flux density at the 
surface of the white dwarf. The relative contribution of the 
two model spectra is determined by the radius ratio of the 
two white dwarfs. The MCMC method maximises the pos¬ 
terior probability, equivalent to minimising y 2 , to find the 
best fit. Each data point is weighted by its uncertainty, with 
no additional weight in favour of either the HST or SDSS 
data. 

The free parameters in our model are the temperatures 
T\ and T 2 , the surface gravities log < 71 , log < 72 , the radius ratio 
R 1 /R 2 , a scale factor s = A-nRf/d 2 to account for the dis¬ 
tance d to SDSS J1257+5428, and the maximum reddening 
along our line of sight E(B — V), which is inc orpo rated us¬ 
ing the expressions presented in lSeatonl (l979l ) an -1 |TT owar thl 
dl983lh We included a uniform prior on the reddening, con¬ 
straining it to 0 < E(B — V)/mag < 0.0173, where the max¬ 
imum is given by the dust map of ISchlaflv fc Finkbeinerl 
(l 201 l|) and we assume a minimum of zero. All other pa¬ 
rameters are left unconstrained. We chose not to include a 
prior on the distance based on the parallax measurements, 
to allow a self-consistency check afterwards. 

The results presented here are based on converged 
chains, from which the so-called burn-in phase is removed. 
We have also thinned the chains, by only storing each 20 th 
model, in order to remove any correlation that may be 
present between subsequent models in the unthinned chain. 


4.1 The hot, massive white dwarf and possible 
pulsations 

With an effective temperature of T 2 = 13030 ± 70 K 
and a surface gravity of log (72 = 8.73 ± 0.05, detailed 
evolutionary models show that the secondary star has 
a mass of M 2 = 1.06 ± 0.05 Mg. The corresponding 
cooling age is T 2 = 1.0 Gyr or 1.2 Gyr, with an esti¬ 
mated uncertainty of 0.1 Gyr, for carbon/ oxygen and oxy¬ 
gen/neon white dwarf models respectively (I Trem blay e t al.l 


l201ll ; iKowalski fc Saumonll2006l : lAlthaus et al. 2007l 1 2 l. The 


values for the mass and surface gravity translate into a ra¬ 
dius of i ?2 = 0.0074 ± 0.0006 Rg. These results are in 
agreement with those of iKulkarni fe van Kerkwiikl J2010ll . 
and the results of the fits to the phase resolved and 
ultraviolet-op tical spectral energy distribution presented in 
iMarsh et al.l (1201 lh . 

At different composition-dependent epochs during a 
white dwarf’s cooling process, the star experiences non- 
radial gravity-mode pulsations. The atmospheric parame¬ 
ters of the secondary white dwarf place it inside the em¬ 
pirical and theoretical instability strip for white dwarfs 
with hydrogen-ric h atm ospheres dGianninas et al] l2014al : 
IVan Grootel et ahl [2013). At this high surface gravity, 
there are only two confirmed white dwarf pulsators 
dGianninas et al . 201ll ; iHermes et al.l l2013l l . Note that this 
empirical instability strip is based on atmospheric param¬ 
eters determined from Balmer line fits, in which the mod¬ 
els used include a ID mixing-length theory to approximate 
convective motion. The most recent models are based on 
3D simulations instead, and give slightly different results 
for both white dwar f temperatures and surface gravities 
dTremblav et al.ll2013l f. Because our atmospheric parameters 
were not obtained through Balmer line fits, these do not suf¬ 
fer from inaccuracies in the ID models. To facilitate direct 
comparison with the empirical instability strip, we there¬ 
fore decided to ‘correct’ our results using the offsets from 


1 http: //www.sdss.org/drl2/algorithms/fluxcal/#SDSStoAB 


2 See http://www.astro.umontreal.ca/~bergeron/CoolingModels 
and http://fcaglp.fcaglp.unlp.edu.ar/evolgroup/. 
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Table 1. White dwarf parameter results from the MCMC analysis performed on the HST+ COS, HST+STIS and SDSS data. The 
reddening is constrained to 0 < E(B — V) < 0.0173 by a uniform prior. Numbers in parentheses indicate statistical uncertainties in the 
last digit(s). The distance is calculated from the scale factor s (see Section [3] for details). The cooling ages t 2 in each column are based 
on carbon/oxygen and oxygen/neon white dwarf models, with an estimated uncertainty of 0.1 Gyr. 


parameters 

MCMC results 

best model 

MCMC 

MCMC 

MCMC 



(see Fig. O 

(fixed log gi ) 

(fixed log g i) 

(fixed log gi ) 

T 2 (K) 

13030(70) 

13033 

13050(59) 

12965(82) 

12811(94) 

Ti (K) 

6400(37) 

6399 

6402(38) 

6395(29) 

6460(24) 

E(B - V) (mag) 

0.0089(34) 

0.0101 

0.0109(23) 

0.0038(21) 

0.0008(8) 

l°g(32) 

8.73(5) 

8.72 

8.73(5) 

8.70(7) 

8.61(9) 

log(si) 

5.26(36) 

5.10 

5.0 

6.0 

7.0 

R 1 /R 2 

4.27(9) 

4.27 

4.29(9) 

4.21(10) 

3.89(9) 

d (pc) 

102(9) 

103 

102(9) 

105(8) 

109(8) 

M 2 (Mq) 

1.06(5) 

1.05 

1.06(5) 

1.04(5) 

1.00(5) 

t 2 (Gyr) 

1.0 / 1.2 

1.0 / 1.2 

1.0 / 1.2 

1.0 / 1.2 

0.9 / 1.2 

Details of the various 

; fits 





minimum 

5771(4) 

5764 

5771(4) 

5774(4) 

5817(4) 

degrees of freedom 

11071 

11071 

11072 

11072 

11072 


3D to ID parameters llTremblav et al]l2013t AT ~ 250 K; 
Alog g ~ 0.01), rather than correcting every other source 
from ID to 3D. 

Although the secondary is placed ~ 600 K from the 
blue edge inside the instability strip, we did not detect any 
pulsations in the time-tagged HST COS data down to an 
amplitude of 1.7%, equivalent to 18 mmag at the 3<r limit. 
To bring this limit down, we obtained the ULTRASPEC 
data. However, these also do not show any pulsations with 
an amplitude exceeding 0.5%. In the g' band light curve 
the contribution of the secondary white dwarf is diluted by 
that of the primary, as the latter contributes 1.6 times as 
much flux at these wavelengths. This puts the 3<r pulsa¬ 
tion amplitude limit at 14 mmag. Note that the HST limit 
is from data at far-ultraviolet wavelengths, where pulsation 
amplitu des are generally muc h larger than at optical wave¬ 
lengths dRobinson et alll995l ), and may therefore still be the 
stronger limit even though the absolute value is somewhat 
higher than that from the ULTRASPEC data. Pulsation am¬ 
plitudes tend to decline for w hite dwarfs with effect ive tem¬ 
peratures exceeding 11500 K (iMukadam et alj|2006l ). and so 
it is possible that they are still present, but with amplitudes 
below the limits presented here. 


4.2 The cool low-mass white dwarf 

The secondary white dwarf mass determined above com¬ 
bined with t he radial velocity v ariation of K\ = 330 km/s 
measured bv lMarsh et al.1 (1201 ll ) put an upper limit on the 
mass of the primary white dwarf at Mi ^ 0.24 Mq (see 
their Fig. 6). This is consistent with the system not being a 
supernova Type la progenitor, as we ll as with the favoured 
solution found in lMarsh et alj d201ll ). 

One interesting result from our analysis is that the data 
strongly suggest that the surface gravity of the primary, 
cooler white dwarf is close to log g = 5.3 (see Fig. |T|). How¬ 
ever, given the radius ratio of R 1 /R 2 = 4.27 and a maxi¬ 
mum possible mass ratio of M 2 /Mi ~ 10, gene rously assum¬ 
ing a minimum white dwa rf mass of 0.1 Mq dlstrate et aid 
l2014al : lAlthaus et aid 120131 ) , the surface gravities can differ 


by log <?2 - log gi ~ 2.3 at most. Given that the surface 
gravity of the hot white dwarf is well-constrained by the 
features in the HST+ COS data we therefore believe that 
the surface gravity of the cool white dwarf should be closer 
to log g\ ~ 6.5. In addition, there is no indication of any ab¬ 
sorption lines besides the Balmer lines, even though at the 
very least the Ca H /K lines are often present in white dwarfs 
with log g < 6 dHermes et aid l20 l ll: iBrown et all 12013 ; 


iGianninas et all l2014bl : iKaplan et al.ll2013l ) . This therefore 

also points towards a surface gravity larger than 5.3 for the 
low-mass white dwarf in SDSS J1257+5428. We do not know 
why the data imply the low surface gravity we find in an 
unconstrained fit. Considering the entire range of possible 
white dwarf surface gravities, a log gi ~ 6.5 is still at the 
low end, and the combination with the low effective temper¬ 
ature is unprecedented, making it difficult to draw robust 
conclusions. 


For these reasons, we reanalysed the data while keeping 
log gi fixed, choosing values of 5.0, 6.0 and 7.0. The results 
are listed in the last three columns of Table [T] The large 
changes in log <?i have relatively little effect on the \ 2 value. 
The main difference between these results and those from 
our initial MCMC is in the values of the reddening and the 
radius ratio. The reddening decreases significantly, becom¬ 
ing consistent with zero when the primary white dwarf’s 
surface gravity is fixed at higher values. This behaviour is 
likely caused by the near-ultraviolet feature in the interstel¬ 
lar extinction curve, which is adjusted to compensate for 
the change in the cool white dwarf’s spectrum, which starts 
contributing to the total flux in this same wavelength range. 
The variation in the other parameters illustrates the extent 
of the systematic uncertainties, which are ~ 150 K for T 2 , 
~ 50 K for Ti, and ~ 0.05 for log g 2 . Note that these un¬ 
certainties are too small to move the secondary out of the 
instability strip. The best fits from the three MCMC runs 
with fixed, different values of log gi are shown in Fig. [3] 
Comparison of these models wit h the Balme r lines in the 
WHT+ISIS spectra presented in lMarsh et al.l (1201 ll ) shows 
that the model with log gi = 6.0 matches the depths of those 
lines best, consistent with our reasoning above. From now 
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Figure 1. Converged MCMC chain projected onto 2-dimensional parameter spaces and showing histograms for the individual free 
parameters in the fits. The contours are at the la, 2 a and 3a levels, and include 68%, 95% and 99.7% of the data respectively. The 
orange squares and vertical lines indicate the best fitting model with x 2 = 5764, as listed in Tabic [T] 


on we therefore assume that log gi ~ 6.0 - 6.5, which agrees 
with the results from the unconstrained MCMC analysis at 
the ~ 3 <t level (see Fig. [lj. 

Given that the surface gravity of the cool white dwarf 
is poorly constrained by the spectra we do not rely on it 
hereafter, and instead use the radius derived for the hot, 
secondary white dwarf in section 14.11 and the radius ratio 
from the MCMC analysis. The latter is constrained by the 
relative flux contributions of the two white dwarfs across the 
spectral energy distribution, and translates into a radius for 
the primary of R\ = 0.032 ± 0.003 Rq. In Fig. [4] we show 
this value and the effective temperature for the cool white 
dwarf, together with e volution ary m odels for white dwarfs 
of different mass from llstrate et al.l (l2014bll . These models 
were obtained for ELM white dwarfs in close binaries with 
neutron stars, but we expect the white dwarf’s formation via 
Roche lobe overflow and detachment to proceed similarly in¬ 


dependent of the nature of the companion, apart from pos¬ 
sible issues of mass-transfer instability. To avoid cluttering 
the figure, we only selected a few of the many models with 
various values of the initial mass of the donor star (the pro¬ 
genitor of the helium white dwarf), the index of magnetic 
b raking, and the m ass of the neutron star companion (see 
llstrate et alj |2014a for further discussion). Our results in¬ 
dicate that the cool white dwarf has a low mass, close to 
0.2 M@, consistent with a low surface gravity. However, the 
models also show that such low-mass white dwarfs take ^ 5 
Gyrs to cool to a temperature of 6400 K, much longer than 
the cooling age derived for the hot white dwarf, which is 
close to 1 Gyr. These values suggest, surprisingly, that the 
low-mass white dwarf formed first. 

Fig. 0 shows a larger area of the same parameter space 
as shown in Fig. [4] now also including E LM white dwarf 
cooling models from lAlthaus et al.l (|2013l l . It is clear from 
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Figure 2. Top panel: best fit model spectra for the double white dwarf binary SDSS J1257+5428 (solid black line) and the individual 
white dwarfs (dotted grey lines). The HST+COS and STIS spectra are shown at 1200 < A/A < 1700 and 1650 < A/A < 3200 
respectively and are binned to 2 A. Solid black dots indicate the SDSS ugriz fluxes (errorbars too small to be seen). The inset shows 
a closer view of the far-UV where the flux is almost entirely dominated by the hot white dwarf. The grey crosses indicate geocoronal 
oxygen emission (A ~ 1300 A), and are excluded from our fits. Bottom panel: residuals of the model SDSS J1257+5428 spectrum folded 
through the SDSS filter curves with respect to the measured SDSS fluxes. 



Figure 3. Top panel: best fit model spectra for the double white dwarf binary SDSS J1257+5428 (solid lines) and the individual white 
dwarfs (dashed and dot-dashed lines) from MCMC fits with log g i fixed at 5.0, 6.0 and 7.0 (black, orange and red respectively). The 
HST+COS and STIS spectra (binned to 2 A in the main panel, and 4 A in the inset) and the SDSS spectrum are shown in grey. Solid 
black dots indicate the SDSS ugriz fluxes (errorbars too small to be seen). The inset highlights the part of the spectrum where the 
models differ most. Bottom panel: residuals of the model SDSS J1257+5428 spectra folded through the SDSS filter curves with respect 
to the measured SDSS fluxes, offset by -100 A, 0, +100 A for log g\ fixed at 7.0, 6.0 and 5.0, respectively. 
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Teff(K) 

Figure 4. Cooling tracks for ELM white dwarfs, together with 
the radius and effective temperature of the cool low-mass white 
dwarf (open square ). The lines are evolutionary models from 
11st rate et al.l (l2014bl see text for more details) for white dwarfs 
of different mass (in Mq, see legend). The various lines are dot¬ 
ted up to a cooling age of 1 Gyr (stars), and solid after. Triangles 
and diamonds are placed at cooling ages of 2.5 and 5 Gyr for each 
track. At T e ff = 6400 K, the white dwarf cooling ages are roughly 
13, 13, 10.5, 7.2 and 6.7 Gyr, with increasing mass, respectively. 



Teff (K) 

Figure 5. Cooling tracks for ELM white dwarfs, together with 
the radius and effective temperature of the cool low-mass white 
dwarf (square). The soli d and dotted lines are evolutionary 
models from l ist rate et al.l d2014bl , labelled I in the legend) and 
lAlthaus et al.l (120131 . labelled A) respectively. For each, a model 
with (labelled F) and without CNO flashes is included. The white 
dwarf masses are in Mq as in the legend. The stars, triangles and 
diamonds indicate cooling ages of 1, 2.5 and 5 Gyr, respectively. 
The horizontal dashed line indicates the size of the ELM white 
dwarf’s Roche lobe in the current binary configuration. 


this figure that the ELM white dwarf in SDSS J1257+5428 
has settled on the cooling track and is not currently in a 
CNO flash cycle. Only ELM white dwarfs that exceed a 
certain mass experience CNO flashes, during which the thick 
hydrogen layer is quickly consumed, thereby speeding up the 
entire cooling process. Our upper limi t of 0.24 Mq is just 
above the minimum mass of 0.1 8 Mq (lAlthaus et al.|[2013l l 
- 0.20 Mq dlstrate et al.ll2014blf necessary for cooling with 
CNO flashes. 

As demonstrated in Figs. [4] and [5] the age of the pri¬ 
mary white dwarf estimated from current cooling models is 
very sensitive to both its mass and the degree of element dif¬ 
fusion. The shortest possible cooli ng age for the E LM white 
dwarf is given by a model from lAlthaus et al.l d2013ll . in 
which the white dwarf is formed with a mass of 0.203 Mq, 
experiences CNO flashes, and takes 1.6 Gyr to reach a tem- 
per ature of 6400 K. The diffe rence in cooling ages between 
the lAlthaus et aid d2013l l and llstrate et aid (l2014bl) models 
are most likely related to the amount of element diffusion 
(for example, via gravita tional settling and radiative levita¬ 
tion, (Althaus_et_alJl200l|). The former models are calculated 
with gravitational settling, whereas the latter models do not 
include this effect. In addition, the treatment of convection 
may play a role. Finally, long-term helium white dwarf cool¬ 
ing (beyond the proto-white dwarf stage) could also be af¬ 
fected by rotation of the white dwarf, which might lead to 
significant mixing and thus prevention of strong element dif¬ 
fusion. New models investigating these issues are currently 
in progress (Istrate et al., in pr ep). _ _ 

The shorter cooling age of lAlthaus et al.l (l2013h is still 
too long to resolve the paradox of the formation of this 
binary. This is compounded by the time the ELM white 
dwarf took to form, since its progenito r most likely had 
a mass < 1.6 Mq dlstrate et al.ll2014a), and thus had a 
main-sequence lifetime of order 1.5 Gyr ( Hurley et alJl2000h . 
which needs to be added to the white dwarf cooling age to es¬ 
timate its total age. Therefore it appears impossible to avoid 


the conclusion that the ELM white dwarf is older than its 
massive white dwarf companion. 


4.3 Distance to SDSS J1257+5428 

Using the results for the scale factor, the radius ratio and 
the secondary radius from our MCMC analysis, we are able 
to calculate the distance via d = Riy/An/s for which we 
find d = 105 ± 8 pc. This is consistent with the distance 
derived using parallax measurements, where d = pc, 

as presented in Sect. 12.21 indicating that our analysis is sen¬ 
sible. At this point we could redo our analysis and include 
a prior on the distance, based on the parallax observations. 
However, given that the uncertainty on our current result 
is smaller than that from the parallax measurements, the 
prior would have little effect. Given furthermore that the 
scale factor s does not correlate significantly with any of the 
other free parameters, the values of these free parameters 
would change little and so we refrain from reanalysing the 
data. 


5 DISCUSSION 


The combined results of the HST data and evolutionary 
models for low-mass white dwarfs present us with an intrigu¬ 
ing puzzle. The secondary white dwarf has a mass just over 
1 Mq, which is near the threshold separating white dwarfs 
with a chemical core-composition of carbo n/oxygcn from 
those with oxygen/neon dominated cores dLazarus et al.l 
I2014T) . If it was an isolated star, we could use an initial- 
final mass relation to obta in an initia l z ero-ag e main se¬ 
quence mass of 5 - 6 Mq dCatalan et alj|2008h. for which 
main -sequence lifetimes are close to 100 Myr ( Hurley et al.l 
l200d) . In close binaries the initial masses are often greater 
than those predicted from initial-final mass relations due to 
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interactions between the stars, and so these numbers have 
to be considered cautiously. Nonetheless, together with the 
cooling age of ~ 1 Gyr, it allows us to estimate the total age 
of the hot, massive white dwarf as 1.1 ± 0.1 Gyr. The low 
temperature of the primary, low-mass white dwarf combined 
with evolutionary models shows that the age of the primary 
white dwarf is at least ~ 1.6 - 5 Gyr. Given the 1.6 M@ 
maximum progenitor mass, its total age is > 3 Gyr. 

We considered whether the cooling age of the massive 
white dwarf could have been reset by accretion heating dur¬ 
ing the formation of the ELM white dwarf companion. How¬ 
ever, this would imply that its cooling age would now be the 
same or longer (if only partially reset) than that of the ELM 
white dwarf, which does not explain what we see. There 
should be accretion after the birth of the ELM white dwarf 
during the CNO flashes as the white dwarf fills its Roche 
lobe (see Fig. E3>. However, these events are very short lived 
(~ 100 yr) and cannot significantly alter the thermal struc¬ 
ture of the massive white d warf which takes ~ 10 6 yr to 
change (|Bildsten et alj|2006h . 

A more exotic possibility is that the massive white 
dwarf formed out of a merger of two white dwarfs roughly 
1 Gyr ago, and 4 Gyr after the formation of the ELM white 
dwarf. The pair had to form well before the ELM white 
dwarf and therefore survive at least 4 Gy r before merg¬ 
ing. Considerations of dynamical stability (lEggleton et al.l 
I1989l l show that if the outer period of this hypothetical 
triple matched today’s 4.6 h period the inner period would 
have had to have been < 1 h. This would result in a merger 
timescale well short of the 4 Gyr minimum. Therefore the 
triple scenario also requires shrinkage of the outermost or¬ 
bit, which implies that the merger was a common-envelope 
event that shrunk both the inner binary and the outer bi¬ 
nary / ELM white dwarf orbit. We cannot say whether this 
is impossible, but it seems unlikely; simulations of white 
dwarf mergers seem to show that the merged object doe s 
not expand significantly (IShen et al .11*201 j ; I Dan et al.ll201llh 
If anything, one might expect that angular momentum from 
the merged pair would be transferred to the outermost orbit, 
resulting in a period increase, not the necessary decrease. 
Even if the proposed scenario is possible, it is hard to see 
how an initial configuration of a tight inner binary contain¬ 
ing at least one carbon/oxygen white dwarf in a close triple 
with an ELM white dwarf could have formed. 

Finally, it is possible that SDSS J1257+5428 is not 
a close double white dwarf, contrary to our assumption 
throughout the analysis presented here. As it has not 
been ruled out that the broad Balmer lines f rom the sec¬ 
onda r y massive white dwarf are sta tionary dMarsh et al.l 
1201 ll ; Kulkarni fe van Kerkwiikll2010l b the system could be 
a triple or the massive white dwarf could be aligned per 
chance with the ELM WD binary instead. Perhaps the low- 
mass white dwarf is in a close binary with an unseen massive 
companion such as a neutron star, while the hotter, massive 
white dwarf is a wide companion. Recently, I Ransom et al.l 
(l2014h discovered a triple system in the Galactic disk con¬ 
sisting of a neutron star and two white dwarfs, of which one 
is very low mass, and hence natur e is apparently producing 
such triple compact star systems dTauris &: van den Heuvell 
120141 1. However, in this scenario the problem with the in¬ 
compatible cooling ages and masses remains, unless the hot 
white dwarf was captured as the third component later on 


and did not form at the same time as the close binary. Such 
an unusual scenario is only likely within a dense stellar clus¬ 
ter environment. Inspection of the HST +STIS acquisition 
image reveals that the point-spread function from the source 
is consistent with being a point source, and so the stars 
would have to be extremely well aligned if it was a chance 
alignment. This is also an argument against the system be¬ 
ing a wide triple, although a close multiple system with a 
separation <10 AU at the time of the observations cannot 
be ruled out. 

Irrespective of the above possibilities, any binary or 
triple system, in which both of the observed white dwarfs 
discussed in this paper were formed, is difficult to reconcile 
with binary stellar evolution. This is mainly due to the fact 
that the progenitor star of the low-ma ss helium white dwar f 
most likely had a mass of 1 - 2 Mg Jlstrate et al.ll2014al h 
and thus a much longer nuclear burning timescale compared 
to that of the 5-6 Mq progenitor of the ~ 1 Mg secondary 
massive white dwarf. 

Future observations to clarify the nature of 
SDSS J1257+5428 could include radio observations to 
search for a neutron star component, as well as phase- 
resolved spectroscopy to measure (or put an upper limit 
on) the radial velocity of the massive white dwarf. If such 
observations confirm the common binary nature of the two 
white dwarfs investigated here, we might be able to use 
their measured masses, radii and temperatures to constrain 
binary evolution and white dwarf cooling models further. 


6 CONCLUSIONS 

We have analysed the spectral energy distribution of the 
double white dwarf SDSS J1257+5428, consisting of HST 
COS and STIS data and ugriz flux measurements from 
SDSS. The effective temperature and surface gravity of the 
hot white dwarf are found to be T 2 = 13030 ± 70 ± 150 K 
and log g 2 = 8.73 ± 0.05 ± 0.05. Evolutionary models show 
that this white dwarf has a mass of M 2 = 1.06 Mg and 
a cooling age of T 2 ~ 1 Gyr. The atmospheric parameters 
place the star inside the ZZ Ceti instability strip, but we did 
not find any pulsations with amplitudes exceeding 18 mmag 
at far-ultraviolet wavelengths or 14 mmag in the optical g' 
band. 

The temperature for the cool white dwarf is T\ = 6400 
± 37 db 50 K, while its surface gravity is constrained to 
log g\ ~ 6.0 - 6.5 by the radius ratio (in turn constrained 
by the relative flux contributions of the two white dwarfs), 
yieldin g a best ma s s esti mate of ^ 0.24 Mg, in agreement 
with iMarsh et al.l (1201 lh . Using evolutionary models we 
find that the age must be > 3 Gyrs, significantly longer 
than the 1.1 Gyr age of the hot white dwarf. The odd com¬ 
bination of both a higher temperature and a higher mass 
for the secondary white dwarf thus cannot be explained by 
substantial accretion during the time the primary white 
dwarf’s progenitor evolved. The difference in cooling ages 
also rules out recent accretion-induced heating as the cause 
of the significant temperature difference between these two 
white dwarfs. Therefore the data surprisingly suggest that 
the low-mass progenitor of the primary white dwarf evolved 
before the high-mass progenitor of the secondary white 
dwarf, thus posing an interesting puzzle regarding their 


© 2015 RAS, MNRAS 000. 019] 































A double white dwarf with a paradoxical origin? 9 


formation scenario. 
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